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ABSTRACT 


A  computer  program  incorporating  the  least  squares  stepwise  multiple  regression 
procedure  has  been  used  to  derive  polynomial  equations  for  converting  coordinates  from 
one  geodetic  datim  to  another  (datum  shifts!  as  a  function  of  latitude  and  longitude. 
Reference  coordinate  differences  throughout  the  geographic  area  or  datum  are  used  for 
the  regression.  Therefore,  the  regression  exhibits  sensitivity  to  regional  variations  in 
the  coordinate  differences.  (These  variations  are  closely  related  to  changes  in  the 
geodetic  control  in  localized  areas.!  The  resulting  equations  reflect  this  in  precision. 
Contrariwise,  the  traditionally  used  Molodensky  Coordinate  Transformation  Formulas, 
as  conventionally  applied  over  a  geographic  area,  utilize  area  average  rectangular 
coordinate  differences  (AX,  AY,  AZL 

The  stepwise  regression  program  evaluates  as  many  variables  as  desired  (typically 
1001;  e.g.,  U7  V5  is  a  single  variable,  where  U  and  V  are  normalized  latitude  and 
longitude,  respectively.  It  sequentially  incorporates  into  the  equation  the  variable 
providing  the  most  improvement  in  fitting  the  reference  coordinate  differences.  After 
such  incorporation  of  a  variable,  all  variables  previously  incorporated  into  the  equation 
are  examined,  and  any  no  longer  significant  are  removed.  This  continues  until  an 
equation  of  specified  precision  is  obtained.  Thus,  an  equation  of  given  precision 
contains  a  minimum  number  of  terms  and  is  therefore  relatively  computer  efficient. 
This  has  particular  relevance  to  near  real  time  applications  involving  computer  storage 
and  computation  time  constraints. 


INTRODUCTION 


During  the  last  two  decades,  considerable  progress  has  been  made  within  the 
Department  of  Defense  (DoD)  in  reducing  the  number  of  geodetic  datums  used  for 
mapping,  charting,  and  geodetic  (MC&G!  purposes.  This  has  been  accomplished  through 
the  development  and  use  of  various  world  geodetic  systems,  the  latest  being  DoD  World 
Geodetic  System  1972  (WGS  72\  (Reference  l!.  Although  the  number  of  geodetic 
datums  actively  used  by  the  DoD  has  decreased  considerably,  multiple  datums  are  still 
in  use  for  many  geographic  areas  and  the  transformation  of  the  geodetic  coordinates  of 
sites  from  one  datum  to  another  is  often  required.  These  coordinate  transformations 
are  accomplished  using  the  Molodensky  Coordinate  Transformation  Formulas.  Another 
technique  for  accomplishing  coordinate  transformations,  multiple  regression  formulas 
(equations),  is  examined  in  this  paper  for  its  applicability  to  provide  precise  coordinate 
transformations  in  near  real  time. 


REGRESSION  PROCEDURE 


The  terms  A<|>,  AX,  AH,  AX,  AY,  and  AZ  represent  reference  coordinate 
differences  (reference  datum  shifts)  as  well  as  regression  equation  computed  datum 
shifts.  This  is  applicable  to  geodetic  coordinates  latitude  (<}>),  longitude  (X),  height  (H), 
and  rectangular  coordinates  X,  Y,  and  Z.  Derived  datum  shift  regression  equations 
provide  datum  shifts  to  the  equation  user  at  any  point  within  the  applicable  regression 
area.  These  equations  are  derived  by  regressing,  for  known  orientation  (reference) 
stations  throughout  the  area,  the  reference  coordinate  difference  dependent  variables 
A<J>,  AX,  AH,  AX,  AY,  and  AZ  individually  on  selective  variables,  sometimes  called 
independent  variables.  For  the  application  subsequently  presented  herein,  these  latter 
variables  consist  of  products  of  powers  of  normalized  latitude  (U)  and  normalized 
longitude  (V);  e.g.,  U3  V1*  is  a  single  variable.  Since  regression  is  on  these  polynomial 
two-dimensional  variables,  the  resulting  regression  equations  contain  these  variables. 

The  basic  idea  of  the  stepwise  multiple  regression  procedure  used  for  this  paper  is 
to  perform  the  regression  as  a  series  of  straight  line  regressions  (steps).  As  many 
variables  as  desired  are  evaluated  in  the  regression,  but  only  a  relatively  small  number 
are  normally  incorporated  in  the  derived  equation.  The  stepwise  multiple  regression 
procedure  sequentially  adds  one  variable  at  a  time  (stepl  to  the  equation;  namely,  the 
variable  that  provides  the  greatest  improvement  in  fitting  the  reference  coordinate 
differences.  After  a  variable  is  added,  all  variables  previously  incorporated  into  the 
equation  are  examined  for  significance,  and  if  any  is  no  longer  significant  it  is  removed 
(another  step).  Correlations  and  F-tests  provide  the  basis  for  entering  and  removing 
variables.  The  correlation  coefficients  are  adjusted  at  each  step.  This  stepwise 
regression  continues  until  statistical  parameters  (F-values),  used  in  the  F-tests,  are 
satisfied;  in  actual  practice,  until  the  desired  equation  precision  is  obtained.  The 
greater  the  number  of  variables  in  the  regression  equation,  the  more  precise  (better 
fitting)  is  the  equation.  Coefficients  for  each  equation  variable,  and  a  constant  term, 
are  determined. 

This  stepwise  addition  and  removal  of  variables  assures  that  only  significant 
variables  are  retained  in  the  final  equation.  Thus,  a  derived  equation  of  given  precision 
contains  a  minimum  number  of  terms  and  therefore  is  relatively  computer  efficient. 
This  has  particular  relevance  to  near  real  time  applications  involving  computer  storage 
and  computation  time  constraints.  For  maximum  computer  efficiency  in  the  compu¬ 
tation  of  datums  shifts,  the  coefficients  and  exponents  of  the  variables  in  the  set(s)  of 
regression  equations  are  utilized  as  data  in  an  efficient  algorithm. 

The  stepwise  multiple  regression  procedure  utilized  herein  enjoys  relatively  high 
favor  in  the  literature.  Alternative  multiple  regression  procedures  include  backward 
elimination,  forward  selection,  stagewise,  all  possible,  and  variations  thereof.  Pro¬ 
cedures  are  described  in  Reference  2. 

REGRESSION  APPLICATION 


Datum  shift  regression  equations  herein  provide  shifts  between  a  local  geodetic 
datum  and  a  world  geodetic  system;  namely,  between  European  Datim  1950  (ED  50)  and 
WGS  72.  The  equations  cover  a  "limited  ED  50  area"  consisting  of  Denmark,  West 
Germany,  Netherlands,  Belgium,  Luxemburg,  and  France,  as  outlined  in  Figure  1. 
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FIGURE  1 


Herein,  latitude  implies  geodetic  latitude,  and  longitude  implies  geodetic  land  geocen¬ 
tric)  longitude.  Table  1  contains  limited  ED  50  area  datum  shift  equations  for  geodetic 
coordinates  latitude  (<j>),  longitude  (X),  and  height  fH),  and  rectangular  coordinates  X,  Y, 
and  Z.  The  equations  were  derived  to  provide  a  maximum  deviation  from  the  reference 
coordinate  differences  of  1.5  meters. 

These  equations  are  used  as  shown  in  Table  1  when  transforming  from  ED  50  to 
WGS  72  (WGS  72  minus  ED  50),  and  are  used  with  a  change  of  sign  when  transforming 
from  WGS  72  to  ED  50  (ED  50  minus  WGS  72'.  The  latter  can  be  readily  accomplished 
by  merely  changing  the  sign  of  the  equation  datum  shift  result.  Subsequently  discussed 
precision  analysis  provides  verification  for  the  dual  direction  utility,  wherein  the 
latitude  and  longitude  on  either  datum  may  be  used  regardless  of  the  direction  of 
transformation.  In  practice,  latitude  and  longitude  values  are  generally  available  on  the 
datum  that  is  being  transformed,  and  are  used. 

Although  variables  composed  of  powers  of  latitude  and  longitude  may  be  used  in 
the  regression,  it  is  generally  more  efficient  and  convenient  to  use  a  normalized 
latitude  (U)  and  a  normalized  longitude  (V)  as  shown  in  Table -1.  For  example,  in  Table  1 
the  normalized  latitude  U  =  3  d>  -  2.61  is  obtained  from  U  =  3  (<J>  -  0.87),  where  0.87  is 
the  approximate  average  latitude  in  radians  over  the  regression  area,  and  3  is  a 
convenient  factor  which  inhibits  large  values  of  equation  coefficients.  Variables 
evaluated  in  the  regressions  total  99  and  consist  of  all  product  combinations  of  U  and  V 
containing  single  digit  exponents,  selections  from  which  were  made  for  the  Table  1 
equations. 

Figure  1,  which  outlines  the  limited  ED  50  area,  also  shows  the  53  observation  points 
at  which  reference  coordinate  differences  were  used  in  the  regressions  to  derive  the 
equations.  Table  2  contains  this  reference  data.  The  first  33  observation  points  are 
Doppler  stations  at  which  reference  differences  consist  of  the  differences  between  ED 
50  coordinates  obtained  from  ground  survey  and  WGS  72  coordinates  obtained  from 
satellite  Doppler  observations.  The  final  20  observation  points  were  located  to  provide 
more  complete  area  coverage,  and  reference  coordinate  differences  thereat  were 
obtained  by  interpolation.  This  interpolation  was  from  the  reference  coordinate 
differences  at  the  33  Doppler  stations  contained  in  Figure  1  and  Table  2  as  well  as 
additional  (not  shown'  Doppler  stations  external  to  the  limited  ED  50  area. 

The  Defense  Mapping  Agency  Aerospace  Center  (DMAAC)  Datum  Transformation 
Stepwise  Regression  Program  was  used  to  derive  the  datum  shift  regression  equations. 
The  reference  Doppler  station  data  was  provided  by  Mr.  J.  F.  Vines,  DMAAC,  who 
obtains  most  of  such  data  from  the  Satellite  Records  Desk  of  the  DMA  Hydrographic 
Topographic  Center,  which  accumulates  it  from  various  agencies  and  governments.  Mr. 
D.  N.  Huber,  DMAAC,  performed  quality  evaluation  of  this  data.  Aforementioned 
interpolated  reference  coordinate  differences  at  20  observation  points  were  inter¬ 
polated  using  datum  shift  contour  charts  prepared  by  Messrs.  Huber  and  Vines  using  a 
Point  Plotting  and  Contouring  Program  developed  by  the  Ohio  State  University.  Mr.  D. 
Holland,  DMAAC,  performed  the  Molodensky  Formula  computations  discussed  at  the 
end  of  this  paper. 

Table  3  is  a  simmary  of  subsequently  discussed  tables  and  contains  deviations  of 
regression  equation  datum  shifts  relative  to  reference  coordinate  differences  for  each 


TABLE  1 

Datum  Shift  Regression  Equations  Between  Limited  ED  50  Area  and  WGS  72 


Limited  ED  50  Area:  Denmark,  West  Germany,  Netherlands,  Belgium,  Luxemburg,  and  France. 
ED 50:  European  Datum  1930.  Identical  to  the  terminology  European  Datum  (ED). 

WGS  72 :  World  Geodetic  System  1972. 


Geodetic  Latitude  Datum  Shift  ( A4)  in  seconds  (") 

2  2  2 

A4(")  =  -  3.17250  +  1.96761  U  +  0.747893  V  -  0.232615  V  ♦  4.68674  U  V 

Geodetic  (and  Geocentric)  Longitude  Datum  Shift  (AX)  in  seconds  (") 

2  2 

AA(")  =  -  5.03S30  -  1.40710  U  +  1.60471  V  -  0.521318  U  +  0.263364  V 

Geodetic  Height  Datum  Shift  ( AH)  in  meters  (m) 

AH(m)  =  47.1915  -  35.1158  U  -  18.2122  V  +  15.8592  U*  +  264.165  U* 

X  Rectangular  Coordinate  Datum  Shift  ( AX)  in  meters  (m) 

AX(m)  =  -  83.7539  -  1.98841  U  ♦  4.26886  V 

Y  Rectangular  Coordinate  Datum  Shift  (AV)  in  meters  (m) 

2 

AY(m)  =  -  107.458  ♦  10.2703  U  -  12.8061  U  V  ♦  46.6228  U  V 
2  Rectangular  Coordinate  Datum  Shift  ( AZ)  in  meters  (m) 

AZ(m)=  -  121.619+  1.72649  V  -  4.42630  V*  -  625.330U*  -  1292.75  U*  v"  +  10532.7U* 


where  U  =  Normalized  latitude  =  3$  -  2.61 

where  $  =  Geodetic  latitude  in  radians. 

V  =  Normalized  longitude  =  3A  -  0.24 

where  X  =  Geodetic  (and  geocentric)  longitude  in  radians. 

Use  positive  X  from  0*  to  ISO*  east  of  Greenwich. 

Use  negative  X  from  0*to  180*  west  of  Greenwich. 

9  and  X  on  either  datum  (ED  50  or  WGS  72)  may  be  used  regardless  of  datum  shift 
direction.  There  is  no  significant  difference  in  the  computed  datum  shift. 

NOTES: 

1.  For  datum  shifts  from  ED  50  to  WGS  72  (WGS  72  minus  ED  50),  use  equations  as  shown 
above. 

2.  For  datum  shifts  from  WGS  72  to  ED  50  (ED  50  minus  WGS  72),  use  above  equations 
with  a  change  of  sign.  This  can  be  readily  accomplished  by  changing  only  the  sign  of  the 
equation  datum  shift  result. 
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TABLE  3 

REGRESSION  EQUATION  DATUM  SHIFT  PRECISION 
Summary  at  53  Observation  Points  Over  Limited  European  Datum  Area 


Coordinate 

Deviation  between  Regression  Equation  Datum  Shifts 
and  Reference  Coordinate  Differences 

*  ED  50 

to  WGS  72 

**  WGS  72  to  ED  50 

RMS 

Meters 

Maximum 

Meters 

RMS 

Meters 

Maximum 

Meters 

Geodetic  Latitude 

0.63 

1.21 

0.63 

1.21 

Geodetic  Longitude 

0.63 

1.38 

0.63 

1.38 

Geodetic  Height 

0.67 

1.34 

0.67 

1.34 

X 

0.53 

1.36 

0.53 

1.36 

Y 

0.60 

1.46 

0.60 

1.46 

Z 

0.57 

1.41 

0.57 

1.41 

Overall 

0.61 

1.46 

0.61 

1.46 

*  ED  50  geodetic  latitude  and  longitude  used  for  regression  equations. 

**  WGS  "72  geodetic  latitude  and  longitude  used  for  regression  equations. 
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coordinate.  Transforming  from  ED  50  to  WG S  72  using  ED  50  geodetic  Jatitudes  and 
longitudes  for  the  equations  results  in  an  overall  root  mean  square  (RMS)  deviation  of 
0.61  meter  and  an  overall  maximum  deviation  of  1.46  meters.  These  same  overall 
deviations  are  obtained  when  transforming  from  WGS  72  to  ED  50  using  WGS  72 
geodetic  latitudes  and  longitudes  for  the  equations. 

For  each  coordinate,  Tables  4  through  9  contain  results  at  each  of  the  53 
observation  points  when  transforming  from  ED  50  to  WGS  72  (WGS  72  minus  ED  50) 
using  ED  50  geodetic  latitudes  and  longitudes  for  the  regression  equations.  The  tables 
contain  the  reference  coordinate  differences,  the  regression  equation  datum  shifts,  and 
the  deviations  of  the  latter  from  the  former.  At  the  bottom  of  each  table  are  listed  the 
RMS  deviation  and  maximum  deviation  for  the  53  points.  The  tables  also  contain 
normalized  weights  used  in  the  equation  derivations. 

Similar  tables,  not  here  presented,  were  obtained  using  the  same  regression 
equations,  with  sign  change,  for  transformation  in  the  opposite  direction;  namely,  from 
WGS  72  to  ED  50  (ED  50  minus  WGS  72)  using  WGS  72  latitudes  and  longitudes  for  the 
equations.  The  maximum  difference  of  results  between  the  two  directions  at  the  53 
points  for  the  six  coordinates  is  0.02  meter.  Thus,  the  same  equations,  with  the 
appropriate  sign,  may  be  used  for  transformation  in  both  directions. 

Table  10  contains  the  Standard  and  Abridged  Molodensky  Coordinate  Trans¬ 
formation  Formulas.  The  Standard  Molodensky  Formulas  are  longer  and  provide  slightly 
better  precision  than  the  Abridged  Molodensky  Formulas.  For  geodetic  coordinates, 
Table  11  contains  a  summary  of  deviations  of  Standard  Molodensky  Formula  datum 
shifts  from  the  previously  discussed  53  reference  coordinate  differences  over  the 
limited  ED  50  area.  A  single  set  of  area  average  rectangular  coordinate  differences 
was  used  as  data  (AX  =  -513.39  meters,  AY  =  -107.27  meters,  and  AZ  =  -122.07  meters). 
Transforming  f  ^m  ED  50  to  WGS  72  results  in  an  overall  RMS  deviation  of  1.58  meters 
and  an  overall  maximum  deviation  of  6.45  meters.  In  comparison,  overall  regression 
equation  deviations  in  Table  3  are  0.61  meter  RMS  and  1.46  meters  maximum.  Thus,  in 
this  illustration,  regression  equation  precision  is  about  three  times  better  than  that  of 
the  Standard  Molodensky  Formulas. 

Since  reference  coordinate  differences  throughout  the  area  are  used  for  the  regres¬ 
sion,  the  regression  exhibits  sensitivity  to  regional  variations  in  the  coordinate  differ¬ 
ences.  (These  variations,  transmitted  from  the  local  ED  50  coordinates,  are  closely 
related  to  changes  in  the  geodesic  control  in  localized  areas.)  The  resulting  datum  shift 
regression  equations  reflect  this  in  precision.  Contrariwise,  the  traditionally  used 
Molodensky  Coordinate  Transformation  Formulas,  as  conventionally  applied  over  a 
geographic  area,  utilize  area  average  rectangular  coordinate  differences  (AX,  AY,  AZ  ). 
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LATITUDE  OATU"  SMI  f  T  ROOT  MEAN  SDUARE  DEVIATION  .020431  SECONDS  .431276  METERS 

LATITUDE  DATUM  SHIfT  MAIIMUM  DEVIATION  .030143  SECONDS  1.21S562  METERS 


LONGITUDE  REFERENCE  VS  REGRESSION  D  RTUN  SHIFTS  FROM  LOCAL  eATUN  TO  VGS  72  (MGS  72  MINUS  LOCAL  DATUM 
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Table  10 

COORDINATE  TRANSFORMATION  FORMULAS 
GEODETIC  DATUM  TO  WGS  72 

A.  The  Stasdaid  M«Ue«A;  rmnht 

A*"  =  {-AX  sin  4  cob  A  -  aY  sin  4  sin  A  +  aZ  cos  a 
+  as  (R^1  sin  4  cos  *)/a 

+  Af  |R„  (a/b)  +  Rjjih  a)jsin  4  cos  .  NR,,  +  H)  sin  I")  1 

AA"  -  I-aX  sin  A  +  aY  cos  a)  •  |(RN  +  H)  cos  4  sin  1"|  "* 

aH  -  AX  cos  4  cos  A  +  aY  cos  4  sin  A  +  AZ  sin  4 

-As  (s/Rn)  +  Af  (b/s)  Rn  sin1 4 

B.  The  Abridge*  MlbM;  Faramlas 

A#"  -  |-AX  sin  4  cos  A  -aY  sin  4  sin  A  +  aZ  cos  4  +  (sAf  +  fAa)  sin  2  4) 

•  (R„  sin  I"|  1 

aa"  -  |-aX  sin  A  +  aY  cos  Af  •  (R,,  cos  *  sin  1"| 

AH  -  AX  cos  4  cos  A  +  aY  cos  4  sin  A  +  AZ  sin  4  +  (sAf  +  fAa)  sin2  *  -  Aa 

f.  Deflaftiea  af  Terns  ia  the  MsMeasfcy  Fsnaslas 

4.  A,  H  -  geodetic  coordinates  (old  ellipsoid) 

4  -  geodetic  latitude.  The  angle  between  the  earth’s  equatorial  plane  and  the  ellip¬ 
soidal  normal  at  a  point  (measured  positive  north  from  the  equator,  negative 

south). 

A  -  geodetic  longitude.  The  angle  between  the  plane  of  the  Greenwich  meridian  and 
the  plane  of  the  geodetic  meridian  of  the  point  (measured  in  the  plane  of  the 
equator,  positive  east  from  Greenwich). 

H  ”  the  distance  of  a  point  from  the  ellipsoid  measured  along  the  ellipsoidal  normal 
through  the  point. 

H  -  N  +*h 

’Indicates  parameters  which  do  not  appear  in  the  Abridged  Molodensky  Formulas. 
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Table  10  (Cont’d) 


N  -  geoid-ellipsoid  separation.  The  distance  of  the  geoid  above 
(  +  N)  or  below  (-N)  the  ellipsoid. 

*h  -  distance  of  a  point  from  the  geoid  (elevation  above  or  below 
mean  sea  level). 

A*.  AX.  aH  -  corrections  to  transform  the  geodetic  coordinates  from  the 
old  datum  to  WGS. 

aX,  aY.  aZ  -  shifts  between  <  lipsoid  centers  of  the  old  datum  and  WGS. 

a  -  semimajor  axis  of  the  old  ellipsoid. 

*b  -  semiminor  axis  of  the  old  ellipsoid. 

*b/a  -  1  -  f 

f  -  flattening  of  the  old  ellipsoid. 

Aa.  Af  -  differences  between  the  parameters  of  the  old  ellipsoid  and 
the  WGS  ellipsoid  (WGS  minus  old). 

e  -  eccentricity. 

e*  -  2f  -  f’ 

Rs  -  radius  of  curvature  in  the  prime  vertical. 

**RS  -  a/fl-e*  sin2  *)'! 

Rm  -  radius  of  curvature  in  the  meridian. 


•*R„  -  ad-cVd-Vam’*)** 


NOTE:  All  A-quantities  are  formed  by  subtracting  old  ellipsoid  values  from  WGS  ellipsoid  values. 


Indicates  parameters  which  do  not  appear  in  the  Abridged  Molodensky  Formulas. 

••For  desk  calculator  computations  involving  commonly  used  ellipsoids,  these  values  are  given  in 
Latitude  Function  Tables,  i.e..  Latitude  Functions  Clarke  1866  Spheroid,  AMSTM  No.  68,  1967. 
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TABLE  11 


STANDARD  MOLODENSKY  COORDINATE  TRANSFORMATION  FORMULA  DATUM  SHIFT  PRECISION 
Summary  at  53  Points  over  Limited  European  Datum  Area 


Coordinate 

Deviation  between  Standard 
Molodensky  Formula  Datum  Shifts 
and  Reference  Coordinate  Differences 

*  ED  50  to  WGS  72 

RMS 

Meters 

Maximum 

Meters 

Geodetic  Latitude 

1.07 

2.81 

Geodetic  Longitude 

2.26 

6.45 

Geodetic  Height 

1.12 

2.77 

Overall 

1.58 

6.45 

*  ED  50  geodetic  latitude  and  longitude,  and  a  single  set  of  area 
average  AX,  aY,  AZ  rectangular  coordinate  differences  used  for 
.  Standard  Molodensky  Formulas. 


